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UWAGI WSTĘPNE 


Funkcja spektralna (widmowa) jest niewątpliwie jedną z najważniej- 
szych charakterystyk słabo stacjonarnych procesów losowych, a wśród 
nich słabo stacjonarnych procesów autoregresyjnych. O szerokim prak- 
tycznym wykorzystaniu wymienionej fukcji decyduje możliwość dekom- 
pozycji na jej podstawie wariancji procesu na składowe odpowiadające 
różnym częstotliwościom. Pomimo iż problemom estymacji funkcji spek- 
tralnej procesów autoregresyjnych poświęcono wiele prac (np. [1], [4], [5]) 
wydaje się, że niektóre z tych problemów nie zostały jednoznacznie roz- 
wiązane. Wymienić tutaj można zagadnienie wyznaczenia takiej długości 
szeregu czasowego, przy której estymatory widma można uznawać za 
wiarygodne ze statystycznego punktu widzenia oraz problem wyboru ro- 
dzaju estymatora w przypadku, gdy nie jest znany rząd procesu auto- 
regresyjnego. W pracy podjęta zostanie próba ustosunkowania się do 
drugiego z wymienionych powyżej problemów. Obok wyników badań 
symulacyjnych omówione zostaną również wybrane metody oceny rzędu 
procesów autoregresyjnych oraz metody estymacji funkcji spektralnej 
tych procesów. 


ESTYMACJA PARAMETRÓW SŁABO STACJONARNYCH 
PROCESÓW AUTOREGRESYJNYCH 


Załóżmy, ze (X; t=0, +1, +2, ...) oznacza słabo stacjonarny i ergo- 
dyczny proces losowy. Proces ten nazywamy procesem autoregresyjnym 
rzędu p, gdy tworzące go zmienne losowe X, czynią zadość następującej 
równości: | 
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X, = oy X, ayXU sd us + appXt_p + et, (1) 
gdzie 
Ve t 0, X1, €2,..) — proces czysto losowy 
o wariancji 82, 
Opi, Ape, -.., Opp — parametry autoregresji. 


Przyjmując, że 
B*X, = Xo k 2 l, 2, es P 


wzór (1) możemy przedstawić w formie 


F(B)X, = s, (2) 
gdzie 
F(B) = 1 — aB — a,B* — ... — appBP. 


Słaba stacjonarność procesu AR(p) oznacza, iż pierwiastki równania 
charakterystycznego F(B) =0 leżą na zewnątrz okręgu o promieniu jed- 
nostkowym. Jak wiadomo przy charakterystyce słabo stacjonarnych pro- 
cesów losowych operuje się często pojęciami funkcji autokowariancji 
i funkcji autokorelacji. Pierwsza z wymienionych wyżej funkcji definio- 
wana jest w następujący sposób: 


y(t) = cov (X, Xa) 0, ET 2 a (3) 
gdzie y(0) oznacza wariancje procesu (X). | 


Funkcja autokorelacji przedstawiana jest w formie: 
elt) =v(1)/y(0), =0, +1, +2,... (4) 


Estymacja parametrów ekonomicznych procesów losowych wymaga przy- 
jęcia założenia o ich ergodyczności. Wynika to oczywiście z faktu, iż 
w przypadku tego typu procesów nie mamy najczęściej możliwości uzys- 
kania ich wielu realizacji. 
Niech (x; t =1, 2, .., n) stanowi n elementową próbę pochodzącą 
z charakteryzującego się słabą stacjonarnością i ergodycznością procesu 
(X; t=0, +1, +2, ..) o zerowej wartości oczekiwanej. W pracy [3] 
wykazano, ze estymator funkcji y(t) charakteryzujący się zgodnością, 
asymptotyczną nieobciążonością oraz najmniejszą wartością błędu średnio- 
kwadratowego wyraża się wzorem: 
n-t 
ci) = n~ XXt, T=0,1,., m (5) 
t=1 


gdzie m (m&n) oznacza punkt odcięcia funkcji autokowariancji. 
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Oceny funkcji autokorelacji wyznaczane są ze wzoru 2 postaci: 
r(t) = ce(9/c(0. T=0,1,..,m (6) 


Obok parametrów y(t) i (t) cennych informacji o przebiegu procesów 
losowych dostarczają również parametry autoregresji. Oceny tych para- 
metrów uzyskiwane są często poprzez rozwiązanie układu równań 
Yule'a-Walkera przedstawianego w następującej formie: 


R(p)a(p) = r(p. poo (7) 
gdzie 
r(0)r(1) ... r(p— 1) 


r(1)r(0) ... r(p— 2) 
R(p) = : 


r(p— 1)r(p —2) ... r(0) 


a (p) s= [apv Sp» ++, App] 7, 
r(p) = [r(1), r(2), ..., r(p)]7. 


W przypadku, gdy macierz R(p) jest nieosobliwa wartości wekto- 
ra a(p) wyznaczane są na podstawie wzoru: 


a(p) = R(p) 'r(p). (8) 


Oceny Yule'a-Walkera wykorzystywane są najczęściej na etapie identy- 
fikacji modeli autoregresyjnych. Znacznie lepszymi statystycznymi włas- 
nościami w porównaniu z ocenami a(p) uzyskanymi na podstawie wzo- 
ru (8) charaktęryzują się oceny wyznaczane w oparciu o metodę Burga *. 
W metodzie tej oceny parametrów ap; i = 1, 2,..,p otrzymywane są 
w wyniku rozwiązania następującego układu równań: 


dS,+2/da,, = 0 (9) 
ap+ı, j = Api — Apts, p+ı Ap, pj+n j=1,2,..,p (10) 
gdzie 
n—p 
Sp+1 = 2(n— p)! [x — apr Xita TT App Xp) t 
t-1 
+ (X x Api Xt—p+1 TOT ozn App x,)*]. 


1 Opis tej metody można znaleźć między innymi w pracy [16]. 
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W pracy [16] zwraca się uwagę na ścisły związek estymatorów Burga 
z estymatorami^ uzyskiwanymi metodą największej wiarygodności °. 
W empirycznej części pracy wykorzystana zostanie metoda Burga. 


AUTOMATYCZNE METODY OCENY RZĘDU 
PROCESÓW AUTOREGRESYJNYCH 


Załóżmy, że parametry autoregresji Qp,, Ap, ..., dpp OSZACOWANO W opar- 
ciu o próbę X,X,..,X, na podstawie metody Burga. Model (1) można 
więc teraz przedstawić w postaci 


X, = ap Xca FaF app Xt-p FE, t=p+1, pk ZN (11) 
gdzie 

Api, Apa, --., App — Oceny parametrów api, Apa, ..., App, 

e — składnik resztowy. 

Jak wiadomo w procesie identyfikacji modeli ekonometrycznych szcze- 
gólne znaczenie przypisywane jest wariancji składnika  resztowego. 
W przypadku modeli autoregresyjnych rzędu p wariancja ta wyznaczana 
jest w oparciu o następujące wzory * 


pis = (n— 2p)* > (Xt — Apr Xi — + — App Xi_p)” (12) 
t=p+1 
lub 
n 
s2 = (n = p $3 (X i pi Miz 7 7 ux 7 App XQ», 
t=p+l 


J. Andel [2] podkreśla, iż ocena rzędu procesu autoregresyjnego dokony- 
wana jedynie na podstawie analizy zachowania funkcji s? (k = 0,1, ..., k) 
nie daje najlepszych rezultatów z tej prostej przyczyny, iż wzrostowi 
rzędu w modelu autoregresyjnym towarzyszy zmniejszanie się z większą 
lub mniejszą regularnością wartości sę (k = 0,1,.., k). Wiadomo jednak, 
że zwiększanie liczby parametrów w modelu autoregresyjnym nie jest 
korzystne, gdyż prowadzi do problemów natury jnterpretacyjnej. Pod- 
stawę tzw. automatycznych metod doboru rzędu procesów autoregresyj- 
nych stanowi funkcja, która z jednej strony uwzględnia wartości sč, zaś 
z drugiej „łagodzi” zmniejszanie się tych wartości wraz ze wzrostem K. 
Wspomniana funkcja przedstawiana jest w następujący sposób: 


2 W pracy [7] podkreśla się, iż metoda największej wiarygodności daje najlepsze 
ze statystycznego punktu widzenia rezultaty. 
3 Zob. [2]. 
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i(k) = s2[1+h(k)], k=0,1,.,k (14) 


gdzie h(k) jest pewną funkcją rosnącą „łagodzącą” zmniejszanie się 
wartości sę wraz ze wzrostem k. | 

Idea oceny rzędu procesu autoregresyjnego polega na wybraniu takiej 
wartości, k, dla której f(k) (0=k=—K) przyjmuje wartość minimalną, czyli 


p=ks=f(k)= min f) (15) 
3 0XISk 
gdzie p oznacza ocenę rzędu p. 


Szeroką prezentację automatycznych metod doboru rzędu procesów 
autoregresyjnych zawarto w pracy [2]. 

W pracy [9], dokonano w oparciu o symulację komputerową porów- 
nania efektywności spotykanych w praktyce metod. Okazało się, iż naj- 
lepsze ze statystycznego punktu widzenia własności należy przypisać 
metodzie Schwarza i metodzie FPE,(k) z a = 4. Podstawę pierwszej 
z wymienionych metod stanowi funkcja o postaci: 


SR(k) = ln s? + k In(n)/n. (16) 


Funkcję kryterium FPE, (k) wprowadzili R. J. Bhanshali i D. Y. Dow- 
nham [6]. Przedstawiana jest ona w następujący sposób: 


FPE, (k) = s? (1 + ak/n). (17) 


Ponieważ obydwie metody charakteryzują się podobną efektywnością, 
do dalszych badań zdecydowano się zastosować tylko jedną z nich, a mia- 
nowicie metodę Schwarza. 


METODY ESTYMACJI FUNKCJI SPEKTRALNEJ 
PROCESÓW AUTOREGRESYJNYCH SKOŃCZONEGO RZĘDU 


Obok funkcji autokowariancji i autokorelacji cenną charakterystyką 
słabo stacjonarnych procesów losowych jest również funkcja spektralna. 
Funkcja ta definiowana jest jako transformata Fouriera funkcji auto- 
kowariancji. W przypadku, gdy transformacie Fouriera poddawana jest 
funkcja autokorelacji otrzymywana jest funkcja gęstości widmowej. Wy- 
mienione powyżej funkcje przedstawiane są w następującej postaci: 


W(t) = > yaje- fEX—1/2, 1/2> (18) 
G- Seem, fE<—1/2, 12> (19) 


gdzie f oznacza częstotliwość. 


382 Adam Góral 


Jak wiadomo praktyczne walory funkcji widmowej wynikają głównie 
z faktu, iż wariancja (X) wyrażana jest w formie: 
1/2 
+(0) = f W(f) dt. (20) 
-1/2 
Rezultatem wzoru (20) jest często powtarzane w literaturze przedmio- .. 
tu stwierdzenie, iż znajomość wartości funkcji spektralnej słabo stacjo- 
narnego procesu losowego umożliwia dekompozycję całkowitej wariancji 
tego procesu na składowe odpowiadające poszczególnym częstotliwościom 
z przedziału <— 1/2, 1/2>. Bezpośrednie zastosowanie wzoru (18) do esty- 
macji funkcji spektralnej prowadzi do estymatora widma o postaci: 
I(f) = y c(rje 120%,  fE<—1/2, 1/2> (21) 


gdzie c(t) oznacza estymator funkcji y(t). 


Estymator I(f) jest estymatorem niezgodnym. Zgodność estymatora 
widma uzyskiwana jest poprzez zastosowanie okien uśredniających oce- 
ny I(f). Wśród wspomnianych okien na uwagę wydają się zasługiwać: 


1. Okno Hanna 


. fi/2[1 + cos(llr/m)] 1=0,1,..,m 
d(t) = k m (22) 
2. Okno Hamminga 
. (0,54 +0,46 cos(llr/m) 1=0,1,..,m 
3. Okno Parzena 
1 — 6r/m*(1 — t/m) 01 [m/2] 
d(t) = 12(1 —/m) [m/2] + 121m (24) 
0 i^m 
4. Okno Bartletta 
_|1— vm 0<T=m 
d(t) = i C= (25) 


Wykorzystywane w praktyce klasyczne estymatory widma przedsta- 
wiane są w postaci: 
w,(f,) = e(0) + 29 c(x) d(x) cos(2ILf1), (26) 
t=1 
gdzie 
wx(f,) — wartość klasycznego estymatora widma w punkcie fj, 
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c(t) — ocena autokowariancji rzędu rt, 
d(t) — okno korelacyjne, 
f, — częstotliwość, dla której wyznaczana jest ocena widma 


(fj = j2m; j=0,1,.., m). 


Statystyczne własności estymatorów w,(f,) odpowiadających poszcze- 
gólnym oknom omówiono m.in. w pracach [12] i [14]. Ciekawe uwagi 
odnośnie do okien uśredniających można znałeźć w [11] oraz [13]. Ważną 
klasę estymatorów funkcji spektralnej stanowią również estymatory 
autoregresyjne. Jeżeli przyjmiemy, że dysponujemy realizacją procesu 
autoregresyjnego rządu p to autoregresyjny estymator widma będziemy 
mogli przedstawić w następujący sposób: 


p 
w,(f) = s2/|1 — DEE -umn| 2 (27) 
I=1 
gdzie: 
s? — wariancja resztowa w modelu autoregresyjnym rzędu p, 
a (l = 1, 2, ..., p) — oceny parametrów autoregresjj w modelu 
AR(p). 


K. N. Berk [5] wykazał zgodność autoregresyjnego estymatora widma 
przy założeniu, ze p*/n — 0 (gdzie n jest liczbą obserwacji w szeregu cza- 
sowym, na podstawie którego szacowane jest widmo). Jedno z twierdzeń 
podanych przez wymienionego powyżej autora głosi, że asymptotyczny 
rozkład autoregresyjnych estymatorów widma jest rozkładem normal- 
nym. 

Rozważając problemy estymacji funkcji spektralnej słabo stacjo- 
narnych procesów losowych warto zwrócić uwagę, iż teoretycy dostar- 
czają nam informacji odnośnie do asymptotycznych własności estymato- 
rów. Długość szeregów czasowych, stanowiących podstawę badań pra- 
ktycznych sprawia jednak, że nie zawsze uwagi teoretyczne można bez- 
pośrednio wykorzystać. Powoduje to, iż pewne problemy z zakresu esty- 
macji widma nie zostały dotychczas jednoznacznie rozstrzygnięte. Prace * 
Parzena i Kromera dowodzą, iż w przypadku, gdy znany jest rząd pro- 
cesu autoregresyjnego, z którego pochodzi dany szereg czasowy autore- 
gresyjne estymatory widma charakteryzują się lepszymi statystycznymi 
własnościami w porównaniu z estymatorami klasycznymi. Wiemy jed- 
nak, że w praktyce najczęściej zmuszeni jesteśmy poddawać ocenie rząd 
procesu, z którego pochodzi dany szereg czasowy. Żadna z metod, które 
do tego celu mogą być wykorzystane nie daje pełnej gwarancji wyboru 
rzeczywistej wartości tego rzędu. Powyższe uwagi wydają się świadczyć 


4 Informację tę przytoczono za pracą [5]. 
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o konieczności podjęcia badań, które dałyby możliwość odpowiedzi na 
pytanie; czy automatyczny dobór rzędu procesu autoregresyjnego może 
spowodować, iż klasyczne estymatory widma okażą się efektywniejsze 
w porównaniu z autoregresyjnymi. 


ANALIZA PORÓWNAWCZA AUTOREGRESYJNYCH I KLASYCZNYCH 
OCEN WIDMA 


W tej części pracy porównana zostanie efektywność estymatorów 
autoregresyjnych widma z estymatorami klasycznymi w przypadku, gdy 
rząd procesu autoregresyjnego jest znany oraz gdy oceniany jest on 
w oparciu o automatyczną metodę Schwarza. Badania przeprowadzono 
na podstawie realizacji wybranych słabo stacjonarnych i ergodycznych 
procesów AR(1) oraz AR(2). Realizacje procesu AR(1) generowano zgod- 
nie z następującą zależnością rekurencyjną: 


X; =£, 

Xt = Ou Xi + Et, t = 2, 3, LIP) n (28) 
gdzie (s; t = 1, 2, .., n) oznacza ciąg liczb pseudolosowych z rozkładu 

N(0,1). 


Analizę przeprowadzono dla ay = 0,3 i ay, = 0,8. Do generowania cią- 
gów liczb reprezentujących proces AR(2) wykorzystano zależności o po- 
staci: 


Xi = £, 

X = £; 

X, = Oy Xt + Ge Xt +e, t = 2, 8, .., n (29) 
gdzie (s; t = 1,2, .., n) oznacza ciąg liczb pesudolosowych z rozkładu 

N(0,1). 


Badania wykonano na podstawie realizacji AR(2) o parametrach 
az = 0,4iaz = 0,45. Liczby pseudolosowe z rozkładu N(0,1) otrzymywano 
w wyniku zastosowania generatora wykorzystującego centralne twierdze- 
nie graniczne. Opis wymienionego generatora zawarty jest w [17, s. 84]. 
Liczby pseudolosowe z rozkładu równomiernego, niezbędne do uzyskania 
liczb z rozkładu normalnego wyznaczono w oparciu o funkcję RANF(X), 
która jest funkcją standardową w systemie CYBER-72. W celu likwida- 
cji efektu" wynikającego z przyjęcia warunków początkowych w spo- 
sób określony we wzorach (28) i (29), w każdej z wygenerowanych reali- 
zacji pomijano 200 pierwszych obserwacji. W pojedynczym eksperymen- 
cie wyróżnić można następujące etapy postępowania: 


5 Zagadnienie to zostało omówione przez G. S. Fishmana [8]. 


O estymacji funkcji spektralnej procesów autoregresyjnych 385 


1) generowanie realizacji danego procesu losowego, 

2) pominiecie 200 pierwszych obserwacji, 

3) wyznaczenie teoretycznych wartości widma analizowanego procesu 
na podstawie wzoru (27) przy założeniu, że s? oraz ay; 1= 1, 2, ..., 
p zastąpiono teoretycznymi wartościami odpowiednich parametrów, 

4) wyznaczenie autoregresyjnych ocen widma poprzez bezpośrednie 
zastosowanie wzoru (27), 

5) obliczenie wartości łącznego błędu średniokwadratowego estyma- 
tora autoregresy jnego, 

6) wyznaczenie ocen widma metodą standardową na podstawie wzo- 
ru (26) i przy wykorzystaniu wszystkich omówionych we wcześniejszych 
rozważaniach okien uśredniających, 

7) obliczenie wartości łącznego błędu średniokwadratowego dla po- 
szczególnych estymatorów klasycznych, 

8) porównanie wartości łącznego błędu średniokwadratowego estyma- 
tora autoregresyjnego z błędami średniokwadratowymi estymatorów kla- 
sycznych, 

9) 100-krotne powtórzenie czynności 1—8, 

10) wyznaczenie wartości średnich i błędów średniokwadratowych 
wszystkich analizowanych estymatorów, dla każdej częstotliwości od- 
dzielnie. 

Niektóre z przedstawionych wyżej punktów wymagają rozwinięcia. 
W przypadku wyznaczania autoregresyjnych ocen widma eksperyment 
przeprowadzono przy założeniu, że rząd procesu jest znany bądź szaco- 
wany na podstawie metody Schwarza. W obydwu przypadkach wartości 
parametrów autoregresji obliczano na podstawie metody Burga. 

Wartości łącznego błędu średniokwadratowego poszczególnych esty- 
matorów uzyskiwano na podstawie wzorów: 


bs, = (m + 1)* $ [w,(f)) — wa(f,)]* (30) 
j=0 
bsk = (m + 1)” Y [wx(£,) — wolf)’, (31) 
j=0 
gdzie: 
(m + 1) — liczba punktów estymacji widma, 


f, = j/2m (j = 0,1, ..., m) — j-ta częstotliwość, 

w-(f;) — wartość widma teoretycznego uzyskiwana na podstawie 
wzoru (27) przy założeniu, iż przyjęte zostaną w nim 
rzeczywiste wartości poszczególnych parametrów, 
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w,(f,) — wartość autoregresyjnego estymatora widma, 
w,(f;) — wartość klasycznego estymatora widma. 
W celu uzyskania większej przejrzystości wyników analizie poddano 
wartości średnie i błędy $redniokwadratowe następujących statystyk: 


za(f,) = wa(f,)/Wz(£;), j = 0,1, . m (32) 
Z,(fj) = wx(f)/ wf), j=01,..,m (33) 
gdzie oznaczenia są analogiczne do przyjętych poprzednio. 


Oceny wymienionych parametrów wyznaczono na podstawie wzorów 
o postaci: 


100 
z,(£,) = 100 ) z,(£) j = 0,1, ..., m (34) 
i-1 
100 
Zy(t;) = 100+" ze)  j-—01.,m (35) 
i=1 
A 100 
mse [z,(£)] = 1007 V" [n(f) — ]5,  j=0,1,..,m (36) 
i=l 
100 
mse [zy(f;)] = w. [n(f) — ]9,  j—91.,m (37) 
gdzie: 
Z4(f,) — wartość statystyki z,(f,) uzyskana w i-tym powtórzeniu 
eksperymentu, 
Zy(fj) — wartość statystyki z,(f,) uzyskana w i-tym powtórzeniu 
eksperymentu. 


Eksperyment przeprowadzono dla różnej liczby punktów estymacji 
(m) i w oparciu o realizacje o różnych długościach (n). 

Rezultaty omówionych badań prezentowane są w tabelach 1, 2, 3 i 4. 

Liczby zawarte w trzech pierwszych tabelach określają, w ilu przy- 
padkach na 100 powtórzeń eksperymentu estymatory autoregresyjne oka- 
zywały się ,lepsze" (ze względu na wartość łącznego błędu średniokwa- 
dratowego) od estymatorów klasycznych. 

Analiza uzyskanych wyników pozwala przecząco odpowiedzieć na 
postawione w pracy pytanie odnośnie do wpływu automatycznego do- 
boru rzędu procesu autoregresyjnego na efektywność autoregresyjnych 
estymatorów widma w stosunku do estymatorów klasycznych. 

Możemy wprawdzie zauważyć, że w wielu badanych przypadkach 
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Tab. 1. Porównanie autoregresyjnych i klasycznych estymatorów widma 
(AR(1) z ay = 0,3) 

A comparison of autoregressive and classical estimators of the spectrum 
(ARQ) with a = 0,3) 


, Liczba Rodzaj Nazwa estymatora klasycznego 
Liczba punktów 
obserwacji | estymacji "A need matora H Ham- Bart- 
widma utoregresyjnego ann mine lett Parzen 
estymator Ia 69 © 70 58 44 
50 11 estymator II» 87 89 85 77 
estymator I 79 81 67 58 
50 13 estymator II 92 94 90 83 
estymator I 91 93 79 77 
50 17 estymator II 96 96 94 94 
estymator I 96 99 94 90 
100 21 estymator II 98 99 96 95 
estymator I 100 100 97 96 
100 26 estymator II 100 100 99 98 
estymator I 100 100 100 100 
100 34 estymator II 100 100 100 100. 
estymator I 95 95 90 85 
200 21 estymator II 99 99 98 98 
estymator I 100 100 100 100 
200 41 estymator II 100 100 100 100 
estymator I 100 100 100 99 
300 31 estymator II 100 100 100 100 
estymator I 100 100 100 100 
300 38 estymator II 100 100 100 100 


a Estymator I oznacza estymator autoregresyjny, dla którego rzad procesu oceniano 
na podstawie metody Schwarza 

b Estymator II oznacza estymator autoregresyjny, który nie wymagat oceny rzedu 
procesu (rząd był znany) 

c Liczba ta oznacza, iż w 69 przypadkach (na 100) estymator autoregresyjny typu 
I okazał się lepszy (w sensie wartości łącznego błędu średniodwadratowego) od 
estymatora klasycznego wykorzystującego ekno Hanna. 

Źródło: obliczenia własne. 


automatyczny dobór rzędu procesu autoregresyjnego powodował, że esty- 
matory klasyczne okazywały się ,lepsze" od autoregresyjnych dla wię- 
kszej liczby eksperymentów, niż to miało miejsce w sytuacji, gdy rząd 
procesu był znany, lecz wyniki uzyskane w obydwu rozważanych przy- 
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Tab. 2. Porównanie autoregresyjnych i klasycznych estymatorów widma 


A comparison of autoregressive and classical estimators of the spectrum 


Liczba 
obserwacji 


50 


50 


50 


100 


100 


+ 100 


200 


200 


300 


300 


a, b, c — oznaczenia analogiczne do podanych w tablicy 1. 
Źródło: obliczenia własne. 


Liczba 
punktów 
estymacji 

widma 


21 


26 


34 


21 


41 


31 


38 


Adam Góral 
p O 2 mie DA ECC 


(AR(1) z a4, = 0,8) 


(AR(1) with ay = 0,8) 


Rodzaj 
estymatora 


autoregresy jnego 


estymator Ia 


estymator II» 


estymator I 
estymator II 


estymator I 
estymator II 


estymator I 
estymator II 


estymator I 
estymator II 


estymator I 
estymator II 


estymator I 
estymator II 


estymator I 
estymator II 


estymator I 
estymator II 


estymator I 
estymator II 


Nazwa estymatora klasycznego 


Hann 


Ham- 


Bartlett 


Parzen 


——— | q—|\|um i] | 


padkach nie róznily sie od siebie istotnie. Wyjatkowe w stosunku do po- 
wyzszych uwag byly wyniki uzyskane na podstawie realizacji procesu 
AR(1) z au = 0,3. Okazało się bowiem, iż dla małych prób (n = 50) au- 
tomatyczny dobór rzędu procesu autoregresyjnego spowodował, że dla 
m = 10 klasyczny estymator widma wykorzystujący okno Parzena był 
»lepszy" od estymatora autoregresyjnego w 56 powtórzeniach ekspery- 
mentu. Sytuację tę można wytłumaczyć, jeżeli uwzględni się wyniki ba- 
dań dotyczących efektywności metody Schwarza *. Warto również pod- 


* Zob. [9]. 
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Tab. 3. Porównanie autoregresyjnych i klasycznych estymatorów widma 
(AR/2/, Az = 0,4; a22 = 0,45) 

A comparison of autoregressive and classical estimators of the spectrum 
(AR/2/, ag, = 0,4; a2 = 0,45) 


Liczba dial Nazwa estymatora klasycznego 
Liczba punktów m > uż 
obserwacji estymacji ALE y Ham- 

widma estymatora Hann mite Bartlett | Parzen 

estymator Ia 83 ¢ 83 85 85 . 
50 11 estymator IIb 85 84 85 87 
estymator I 81 78 83 85 
50 13 estymator II 83 80 84 87 
estymator I 75 75 80 81 
50 17 estymator II 77 78 81 84 
estymator I 17 73 83 85 
100 21 estymator II 82 81 86 88 
estymator I 72 73 79 79 
100 26 estymator II 78 78 85 85 
estymator I 72 73 77 73 
100 34 estymator II 78 78 79 81 
estymator I 80 80 80 84 
200 21 estymator II 81 81 81 82 
estymator I 76 76 76 77 
200 41 estymator II 75 77 78 78 
estymator I 81 77 81 85 
300 31 estymator II 80 79 84 84 
estymator I 72 69 78 82 
300 38 estymator II 75 74 82 82 


a, b, c — oznaczenia analogiczne do podanych w tablicy 1. 


Źródło: obliczenia własne. 


kreślić, iż w omawianym przypadku dla prób o liczebności n = 100, 200, 
300 autoregresyjne estymatory widma okazywały się lepsze od estyma- 
torów klasycznych we wszystkich 100 powtórzeniach eksperymentu za- 
równo w sytuacji, gdy rząd był znany, jak i w sytuacji, gdy był on oce- 
niany na podstawie metody automatycznej. Zawartość wszystkich czte- 
rech tabel wyraźnie wskazuje, iż estymatory autoregresyjne charaktery- 
zują się lepszymi statystycznymi własnościami w porównaniu z estyma- 
torami klasycznymi. Należy zaznaczyć, że autoregresyjne estymatory 
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Tab. 4. Wartości średnie i błędy $redniokwadratowe estymatorów widma wybran 
The mean values and mean square errors of the estimators of the spectrum of the 


Proces AR (1) 


04, = 0,3 

fi 
z | mse? | z5 | MS€ap zi | msey Za m$ea, 
0 1,024 0,089 1,020 0,054 0,953 0,089 | 1,006 0,186 
1/40 1,017 0,078 0,018 0,051 0,981 0,070 0,961 0,071 
2/40 1,002 0,054 1,012 0,043 1,022 0,067 0,974 0,034 
3/40 0,991 0,035 1,004 0,034 1,030 0,071 1,000 0,034 
4/40 0,990 0,026 0,997 0,026 1,012 0,064 1,015 0,039 
5/40 0,998 0,030 0,992 0,020 0,982 0,056 1,020 0,038 
6/40 1,005 0,037 0,989 0,015 0,964 0,049 1,018 0,032 
7/40 1,004 0,036 0,987 0,012 0,975 0,042 1,013 0,026 
8/40 1,001 0,034 0,986 0,011 1,004 0,042 | 1,007 0,020 
9/40 0,995 0,031 0,986 0,011 1,017 0,053 1,001 0,015 
10/40 0,990 0,027 0,987 0,010 1,005 0,062 0,996 0,013 
11/40 0,985 0,023 0,988 0,011 1,000 0,057 0,993 0,011 
12/40 0,982 0,020 0,990 0,012 1,002 0,051 0,990 0,012 
13/40 0,980 0,018 0,991 0,013 1,000 0,051 0,988 0,013 
14/40 0,981 0,017 - 0,992 0,014 1,008 0,065 0,987 ' 0,014 
15/40 0,983 0,018 0,993 0,014 1,006 0,061 0,987 0,017 
16/40 0,988 0,020 0,994 0,015 0,990 0,043 0,987 0,019 
17/40 0,994 0,025 0,995 0,016 0,993 0,047 0,988 0,022 
18/40 0,999 0,029 0,995 0,016 1,002 0,045 0,989 0,024 
19/40 1,003 0,033 , 0,996 0,016 0,987 0,056 0,989 0,025 


20/40 1,004 0,034 0,996 0,016 0,973 0,085 0,989 0,026 


a — średnia wartość statystyki za (fj przy założeniu, ze estymator autoregresyjny 
b — mse oznacza wartość średniokwadratową odpowiedniego estymatora, 

c — średnia wartość estymatora autoregresyjnego przy założeniu, że rząd procesu 
d — średnia wartość wybranego (najlepszego w danym przypadku) estymatora kl 
Źródło: obliczenia własne. 


widma charakteryzowały się niską obciążonością i niewielkimi wartościa- 
mi błędu średniokwadratowego nawet w przypadku realizacji o 50 ele- 
mentach (ze względu na ograniczoność miejsca wyniki te nie są prezen- 
towane). 


ZAKOŃCZENIE 


Omówione w pracy badania symulacyjne upoważniają do sformuło- 
wania następujących wniosków: 
1) autoregresyjne estymatory widma charakteryzują się „lepszymi” 
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ych procesów autoregresyjnych (n — 200, m — 20) 
autoregressive processes (n — 200, m — 20) 


Proces AR (1) vods AR (2) 


au = 0,8 | = 0,45 

: x ES EXPE ES T : HHE 
Za2. | MSEa2 z mseak (GSA RUE RA z msea MSE€ ae Zk mse, 
1,013 | 0,172. | 0,766 0,148 0,925 0,318 0,940 0,297 0,537 0,276 
0,973 | 0,068 ; 0,969 0,090 0,934 0,057 0,926 0,049 1,230 0,257 
0,969 | 0,025 1,192 0,157 1,004 0,036 0,980 0,023 1,488 0,416 
0,978 | 0,014 1,155 0,148 1,029 0,036 1,003 0,022 1,193 0,177 
0,986 | 0,011 | 1,102 0,129 1,034 0,026 1,013 0,023 1,127 0,146 
0,991 | 0,011 1,028 0,089 1,029 0,033 1,018 0,024 1,023 0,089 
0,994 | 0,010 | 0,978 0,071 1,021 0,029 1,019 0,024 0.977 | 0,067 
0,996 | 0,010 0,977 0,059 1,012 0,026 | 1,019 0,023 0,969 0,060 
0,998 | 06,010 | 1,016 0,062 1004 | 0,024 | 1,018 0,022 1,015 0,065 
0,999 | 0,011 | 1,032 0,077 0,998 0,024 ; 1,016 0,020 1,026 0,078 
1,000 | 0,011 | 0,997 0,091 0,994 0,026 1,013 0,018 0,995 0,090 
1,000 | 0,011 | 0,996 0,077 0,991 0,028 1,009 0,016 1,001 0,075 
1,001 | 0,011 | 1,004 | 0,066 0,990 0,030 | 1,005 0,013 1,009 0,068 

1,001 | 0,011 0,992 0,069 0,991 0,033 0,999 0,011 0,995 0,068 
1,001 | 0,011 1,016 0,096 0,990 0,037 0,993 0,010 1,021 0,100 

1,001 | 0,011 | 1,013 0,085 0,987 0,040 0,987 0,011 1,010 0,081 
1,001 | 0,011 | 0,973 0,059 0,977 0,032 0,982 0,014 0,980 0,059 
1,002 | 0,011 | 0,987 0,068 | 0,968 | 0,024 | 0,980 0,021 0,995 | 0,069 

1,002 | 0,011 | 1,013 0,065 0,975 0,034 0,982 0,032 1,006 0,062 
1,002 | 0,011 | 0,984 | 0,067 | 0,999 | 0,062 | 0,987 | 0,045 | 0,963 | 0,066 
1,002 | 0,011 | 0,952 | 0127 | 1,016 EL Er Roos dod ud aedi | 0,086 0,990 | | m | om | 0,052 | 0,926 | 0,121 


wymaga szacowania rzędu procesu 


jest znany, 
asycznego, 


statystycznymi własnościami w porównaniu z estymatorami klasycznymi 
bez względu na to, czy rząd procesu jest znany, czy też oceniamy go na 
podstawie efektywnej metody automatycznej, l 

2) estymatory autoregresyjne charakteryzują się niską obciążonością 
i wysoką efektywnością nawet w przypadku, gdy ich wartości wyznacza- 
ne są na podstawie szeregów czasowych o długości 50 elementów. 

Na zakończenie należy dodać, iż badania przedstawione w niniejszej 
pracy przeprowadzone zostały w oparciu o programy komputerowe na- 
pisane przez autora. 
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17. 
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PE3OME 


IIpeAtpKHATA NONBITKA IIpoBepKHM rHIIOTe3PI, T.IOCALIEŃ: ABTOMATHAECKUŃ NHOĄÓ0P 


pana aBTOperpeccMOHHOTO mpouecca Tak BIMAET Hà 3HAHEHMA ABTODETPECCHOHHBIX 
OHeHOK CIeKTDa, YTO OHM CTAHOBATCA ,xyZXXe" (B CMbICJIe 3HATEHMA CpenHel KBa- 
nparMuecKoj( OLIKÓKM) Kkuaccwaüeckux oneuok. Bepudukanus STOM runore3bi Obra 
IpoBeqeHa Ha OCHOBe MOJISJMpyiOHIMX MCCJIeJOBaHMÜÁ, OIIMpascb IIa peaJHMSalMio 
cxa6ocranuonapubix M 3proJyMueckgzx mnponeecoB AR (!) u AR (2). PaccMoTpeH1 Tar- 


ace 


upexH aBTOMAaTMHEeCKMX  MeTOJIOB rion6opa Pana aBTODerpecCcHMOHHOTO riporecca 


M nponenypa ONEHKM napaMeTpoB aBToperpeccuu, npegJroxeuHas Byprow. 
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SUMMARY 


The work undertook an attempt to verify the hypothesis that automatic choice 
of the order of the autoregressive process has such an influence on the autoregres- 
sive values of the estimators of the spectrum that, they become "worse" (in the sense 
of the value of mean-square error) that classical estimators. Verification of 
the above hypothesis was made on the basis of simulation investigations, on the 
basis of the realization of weakly stationary, and ergodic processes AR (1) and 
AR (2). The work also discussed the idea of automatic methods of choosing the 
order of the autoregressive process and the procedure of estimating the parameters 
of autoregression as suggested by Burg. 
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